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Abstract 

Mixed spin-1/2 and spin-1 Ising ferrimagnets on a triangular lattice with sublattices A, B and 
C are studied for two spin value distributions (5a, 5b, 5c) = (1/2,1/2,1) and (1/2,1,1) by Monte 
Carlo simulations. The non-bipartite character of the lattice induces geometrical frustration in 
both systems, which leads to the critical behavior rather different from their ferromagnetic coun¬ 
terparts. We confirm second-order phase transitions belonging to the standard Ising universality 
class occurring at higher temperatures, however, in both models these change at tricritical points 
(TCP) to first-order transitions at lower temperatures. In the model (1/2,1/2,1), TCP occurs on 
the boundary between paramagnetic and ferrimagnetic (±1/2, ±1/2, ^1) phases. The boundary 
between two ferrimagnetic phases (±1/2, ±1/2, ^1) and (±1/2, ^1/2,0) at lower temperatures is 
always first order and it is joined by a line of second-order phase transitions between the paramag¬ 
netic and the ferrimagnetic (±1/2, ^1/2, 0) phases at a critical endpoint. The tricritical behavior 
is also confirmed in the model (1/2,1,1) on the boundary between the paramagnetic and ferrimag¬ 
netic (0,±1,±1) phases. 

PACS numbers: 05.50.±q, 64.60.De, 75.10.Hk, 75.30.Kz, 75.50.Gg 

Keywords: Mixed-spin system, Frustrated Ising ferrimagnet, Triangular lattice, Monte Carlo simulation, 
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I. INTRODUCTION 


Mixed-spin Ising systems have been mostly investigated as possible models of some 
types of ferrimagnetic and molecular-based magnetic materials. The used approaches include 
an exact treatment in special cases^^-, mean-held approximatiort^'^, effective-held theory with 
correlations^^—, Monte Carlo simulations^^— and some other methods^^i^— . The main focus 
were their phase diagrams as well as technologically interesting compensation behavior with 
possibility to achieve zero total magnetization by tuning of temperature below the critical 
point. Most of the studies considered the simplest models consisting of two sublattices one 
of which is occupied with spins S' = 1/2 and the other with S' = 1. Such a mixed-spin model 
can be described by the Hamiltonian 

H = ( 1 ) 

{*.i> i 

where cTj = ±1/2 and Sj = ±1,0 are spins on diherent sublattices, {i,j) denotes the sum 
over nearest neighbors, J < 0 is a antiferromagnetic exchange interaction parameter and D 
is a single-ion anisotropy parameter. Negative values of the parameter D favor nonmagnetic 
states with Sj = 0 and positive values magnetic states with Sj = ±1. 

Due to persisting ambiguities majority of the investigations focused on the simplest lat¬ 
tices, i.e., the square in two and cubic in three dimensions. We note that a long standing 
controversy regarding the critical and compensation behaviors even for the most studied 
case of the model on a square lattice was solved only recently by Monte Carlo simulation 
that has convincingly shown^^ that there are neither tricritical nor compensation points, as 
had been suggested by some previous approximative approaches^'^'^'^. On the other hand, 
in the same study the presence of both the tricritical point and a line of compensation 
points was conhrmed in the three-dimensional model on a simple cubic lattice. This hnding 
might suggest that the increased dimensionality is responsible for the appearance of the tri¬ 
critical and compensation behaviors. Nevertheless, our recent study on a triangular lattice 
ferromagnet^ demonstrated that the tricritical point can also appear in a two-dimensional 
lattice as long as the coordination number is sufficiently high. The effect of the coordination 
number in the presence of bond disorder on tricritical behavior of a two dimensional system 
was also recently studied in a random Blume-Capel model on a triangular lattice^. 

We point out that the previous studies were performed on bipartite lattices, in which 
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FIG. 1: (Color online) Mixed-spin S = {Sa, Sb, Sc) models on a triangular lattice consisting of 
sublattices A, B and C, with (a) S = (1/2,1/2,1) mixing - model I and (b) S = (1/2,1,1) mixing 
- model II. Small and large circles denote spin-1/2 and spin-1 sites, respectively. 


case the sign of the exchange interaction is irrelevant to the thermodynamic and critical 
properties of the model in the absence of an external held. On the other hand, the present 
mixed-spin model is considered on a non-bipartite triangnlar lattice, in which case the sign 
of the exchange interaction matters. Namely, in contrast to the ferromagnetic case, the 
ferrimagnetic interaction will indnce geometrical frustration, which can be expected to have 
some impact on the critical behavior. As shown in Fig. [H the lattice consists of three 
sublattices A, B and C, occupied with spins S = (S'a, 5'b, Sc)- This allows to further study 
the model in two different mixing modes. We can consider a mixed-spin S = (1/2,1/2,1) 


model I, as schematically depicted in Fig. 1(a), in which one sublattice is occupied with 
spin S = 1 sites and the remaining two sublattices with spin S' = 1/2 sites. Thus, each 
spin-1 site is surrounded hy z = 6 nearest neighbors with spin S = 1/2. The other way 


of the spin-mixing is realized in a S = (1/2,1,1) model II, shown in Fig. 1(b), which is 
obtained when the spin-1/2 and spin-1 sites in the model I are swapped. The two models 
were shown to display qualitatively different critical behaviors even for the ferromagnetic 
exchange interactions^^. 

The goal of the present study is to examine effects of the geometrical frustration on the 
critical behavior of the above dehned ferrimagnetic mixed-spin systems, to determine their 
phase diagrams and to confront them with their ferromagnetic counterparts as well as the 
pure spin-1/2 and spin-1 antiferromagnetic systems. 
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FIG. 2; (Color online) Time evolutions of the internal energy per site e/| J| and the staggered mag¬ 
netizations msi,ms 3 (see definitions below), starting from the disordered phase at the temperatures 
kBT/\J\ = 0.7 in model I and 1.0 in model II, for D/\J\ = 0 and L = 120. 

II. MONTE CARLO SIMULATION 

In order to study the behavior of various thermodynamic quantities in the parameter 
space and to determine the phase diagrams we use Monte Carlo (MC) simulations with the 
Metropolis update rule and employ the periodic boundary conditions. We consider lattices 
with the size L x L, with L ranging from 24 up to 120. We perform N = 2 x 10^ up 
to 10® MCS (Monte Carlo sweeps), the first 20% of which are used to bring the system 
to equilibrium and then discarded, and the remaining data are used to estimate thermal 
averages and statistical errors. In order to demonstrate that the used MCS is sufficient to 
ensure equilibrium conditions, in Fig. |2] we present MC time evolutions of some relevant 
quantities, such as the order parameters rrisi (model I) and rriss (model II) and the internal 
energy per site e/\J\, starting from the disordered phase. The simulations are performed 
in both models for the largest of the considered lattice sizes L = 120, which is the most 
difficult to equilibrate, at the temperatures in the vicinity of the respective critical points 
for D/\J\ = 0. The plots show that in these cases the equilibrium is reached in less than 
10^ MCS. 

The phase boundaries are roughly determined from the maxima of some thermodynamic 
functions, such as the specihc heat, for a selected hxed value of L. We chose L = 48, as 
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a compromise value above which the specihc heat maxima positions do not change con¬ 
siderably and the phase diagrams can be determined in a relatively wide parameter space 
in a reasonable computational time. In the region where the critical line as a function of 
the single-ion anisotropy parameter D is more or less horizontal it is convenient to obtain 
temperature dependencies of the calculated quantities at a hxed value of D. In such a 
case simulations start from the paramagnetic phase using random initial conhgurations with 
the temperature gradually decreased and a new simulation starting from the hnal conhg- 
uration obtained at the previous temperature. On the other hand, if the phase boundary 
shape changes to vertical we obtain variations of the quantities as functions of the single-ion 
anisotropy parameter D at a hxed temperature. Then simulations start from appropriately 
chosen states (i.e., not necessarily random), expected in the considered region of the param¬ 
eter space. Following the above described approach we ensure that the system is maintained 
close to the equilibrium in the entire range of the changing parameter and thus considerably 
shortens thermalization periods. In order to estimate statistical errors, we perform three 
independent simulations at all considered parameter values. 

At some selected points of the phase boundaries we perform a more thorough hnite- 
size scaling (FSS) analysis in order to determine more precisely the location of the critical 
points and the corresponding critical exponents. In such a case we perform more extensive 
simulations using up to = 10^ MCS and apply the reweighing techniques^. For more 
reliable estimation of statistical errors, in this case we used the F-method^. Having obtained 
the maxima of the relevant quantities, we apply the linear htting procedure for logarithms of 
data with errors, following the method in York et ah—. In order to assess the quality of the 
htting, as a measure of goodness of ht we evaluated an adjusted coefficient of determination^^ 
of the linear ht R^. The critical points and the exponents are then extracted from the FSS 
analysis, using the linear sizes L = 24,48, 72, 96 and 120. 

We calculate the following quantities: the internal energy per spin e = {H)/L‘^, the 
respective sublattice magnetizations per site mx, (X = A, B or C), as order parameters on 
the respective sublattices, which for the model I are given by 


"iA(B) = 3 (|Ma(b)|)/A^ = 3{| CTj 

(2) 

iGA(B) 


mc = 3{\Mc\)/L^ = 3(jj2^^ 

(3) 


iec 
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and for the model II by 


rriA = S{\Ma\)/L^ = 3 (| ^ 4 ^ 

ieA 

mB(c) = 3(|MB(c)|)/i:' = 3(| S,\)/L\ (5) 

ieB(C) 

where (■ ■ ■) denotes thermal average. Based on the ground-state considerations (see below), 
for the identihed ordered phases we additionally dehne the following order parameters for the 
entire system, which take values between 0 in the fully disordered and 1 in the fully ordered 
phase. For the model I we introduce two order parameters (staggered magnetizations per 
site) nisi and 771*2 given by 


m*i = (IM*iI)12 a, + 2 a, - I)/L^, 

ieA jeB fcec 


and 

777*2 = (|M*2|)/F^ = 3(11 '^(Ji 

ieA jeB 

For the model II we dehne the order parameter 777*3 as 


( 6 ) 

(7) 


ms3 = (|M*3|)/L2 = 3(1 (8) 

jeB kec 


Unlike in ferrimagnetic systems on bipartite lattices, the three sublattices in the present 
system facilitate spin arrangements in such a way that the total net magnetization is always 
zero and, therefore, of no practical use. 

Further, we calculate the susceptibilities pertaining to the respective order parameters 
O = Mx (X = A, B, C) and also O = Mgi (7 = 1,2 or 3) 


Xo 


(Q^) - {O? 

NoksT ’ 


(9) 


the specihc heat per site c 


{H^) - {Hy 


NoksT^ ’ ^ ^ 

where No is the number of sites on the (sub)lattice on which O is dehned. Further, we 
dehne the logarithmic derivatives of (O) and (O^) with respect to /3 = l/ksT, 
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( 12 ) 


and finally the fonrth-order Binder cnmnlant Uq corresponding to the qnantity O 


Uo = l 


3 ( 02)5 


(13) 


The above standard thermodynamic qnantities (lUlfTH. as well as the less traditional ones 
dehned by Eqs. (nnm, serve to obtain estimates of the respective critical exponents by FSS 
analysis^*^. In particnlar, we nse the following scaling relations, applied to the maximnm 
valnes of the following fnnctions: 


XO,ma.{L) CX 

(14) 

Cmax{L) OC , 

(15) 

Dlo,max{L) OC , 

(16) 

D2o,max{L) OC , 

(17) 


where a is the critical exponent of the specihc heat and uq-, 7 o are the critical exponents 
of the correlation length and snsceptibility, respectively, pertaining to the qnantity O. In 
case of a hrst-order phase transition, the above qnantities (1141117^ are expected to scale as 
oc L'^, where d = 2 is the system dimension. The order parameter cnmnlant, dehned by 
Eq. f[T5D . can also serve for a simple yet relatively precise location of the phase transition 
point as a point at which the cnmnlant cnrves obtained for different system sizes intersect 
at a nniversal valne, e.g., Uo{Tc) = 0.611 for a two-dimensional Ising model^. 


III. RESULTS 

A. Ground state 

Let ns hrst identify all the possible gronnd states (GS) for entire range of the single-ion 
anisotropy parameter D. Considering the lattice system consisting of three interpenetrating 
snblattices A, B and C, as schematically depicted in Fig.[Tl the Hamiltonians of the respective 
models I and II can be dehned as 

Hi = -j(^ ^ ajSkj-D'^Sl, (18) 

ieAJeB i^A,keC jeB,kGC keC 
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TABLE 1: Ground state configurations and the respective energies for different ranges of the single- 


ion anisotropy parameter. 

Model 

I 

II 

D/\J\ 

State Energy e/| J 

State Energy e/| J| 

(- 00 ,-3/2) FR^ 2 : 

(±1/2, ^1/2,0) -1/4 

P: (0,0,0) 0 


(-3/2, cx)) FR^: (±1/2, ±1/2, :f 1) -3/4-L»/3|J| FR^^: (0,±l,:i:l) -1 - 2L»/3| J| 

E E E sa)-- c'(Es|+ES‘)' 

ieA,j£B ieA,k&C jeB,keC jeB keC 

Focusing on a triangular elementary unit cell consisting of the spins S'a, Sb, Sq, from 
the Hamiltonians flT8|) and flT^ one can obtain expressions for the reduced energies per 
spin e/\J\ of different spin arrangements as functions of D/\J\. Then the ground states 
are determined as conhgurations corresponding to the lowest energies for different values of 
D/\J\, as tabulated in Table [B There are two long-range order (LRO) ferrimagnetic (FR) 
states FR^ and FR 2 in the model I and one LRO ferrimagnetic state FR^^ and one disordered 
paramagnetic (P) phase in the model II. We note that while in FR 2 zero means nonmagnetic 
states {Sj = 0) of spins on sublattice C, in FR^^ zero means magnetic states (a* = ±1/2) of 
spins on sublattice A but equally in states ±1/2 and —1/2, thus giving zero net sublattice 
magnetization. The critical value of the single-ion anisotropy parameter separating the 
respective phases is the same for both models Dc/\J\ = —3/2. 


B. Monte Carlo 


1. Model I: S = {1/2,1/2,1) 


Phase boundaries between the disordered paramagnetic and the respective ordered ferri¬ 
magnetic phases are determined from the specihc heat maxima for a hxed L = 48^ and the 
nature of the ordered phases is established from the introduced order parameters rusi and 


ms 2 - These are shown in Figs. 3(a) and 3(b) for selected values of the parameter D/\J\ below. 


close to and above the critical value Dc/\J\. All the specihc heat curves show pronounced 
sharp peaks, signifying phase transitions to the low-temperature ferrimagnetic states. How¬ 
ever, the (pseudo)transition temperatures appear to be a nonmonotonic functions of D/\J\. 
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FIG. 3: (Color online) Temperature dependencies of (a) the specific heat and (b) the order param¬ 
eters, for selected values of D/\J\ and a fixed lattice size L = 48. In (a) additionally the results 
for L = 72 are shown (squares). 


More specifically, the transition temperatnre for D/\J\ = —1.48, i.e., close to the critical 
value Dc/\J\, is lower than for the other two values below and above Dc/\J\- Moreover, the 
corresponding specific heat maximum has a spike-like shape and its magnitude is one order 
higher then the other maxima, which is typical for a hrst-order phase transition. In order 
to support our claim, that the peaks’ positions do not significantly change above L = 48, 


in Fig. 3(a) we also included the results obtained for L = 72. This will also become evi¬ 
dent later on in the respective phase diagrams in which for some selected points the phase 
transition temperatures estimated for L = 48 will be compared with those determined for L 


extrapolated to inhnity. The order parameters depicted in Fig. 3(b) demonstrate that the 
transition for D/\J\ = —2 is to the state (±1/2, =pl/2, 0), characterized by a finite values 
of ms 2 i while for D/\J\ = —1.48 and 0 the system tends to the state (±1/2, ±1/2, 4=1), 
characterized by a finite values of m^i. The discontinuous behavior of the order parameter 
for D/\J\ = —1.48 corroborates the first-order nature of the transition. 

By FSS analysis at two representative values of D/\J\ = 0 and —2, selected on either 
side of the critical value Dc/\J\ = —3/2, we confirmed that the disorder-to-order phase 
transitions to both (±1/2, ±1/2, =pl) and (±1/2, =|=l/2, 0) ferrimagnetic phases are indeed 
second order and belong to the standard Ising universality class. The slopes of the fitted 
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FIG. 4: (Color online) FSS analysis of the critical exponents ratios l/vo and (a,b), the 

fourth-order cumulant Uq temperature dependencies for different L (c,d) and the Uq data collapse 
analysis (e,f), for O = Mgi at D/\J\ = 0 (left column) and for O = Ms 2 at D/\J\ = —2 (right 
column). The coefficients of determination B? for the respective fits (from top to bottom) are 
0.9999, 0.9998, 0.9997 in (a) and 0.9999, 1.0000, 1.0000 in (b). 
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FIG. 5; (Color online) Order parameter as a fnnction of the increasing (>) and decreasing (<i) 
single-ion anisotropy parameter D/\J\ at varions temperatnres and L = 48. The donble-headed 
arrows mark the hysteresis widths. 


curves in Figs. 4(a) and 4(b) represent ratios of the critical exponents l/i^o and 7o/t'o, 


following from the scaling relations flT^ - (lT7)) . where O = M^i {i = 1,2). We also checked 
that for both D/\J\ = 0 and —2 the specihc heat maxima follow the logarithmic scaling 
Cmax = Co + Ciln(L), as expected for the Ising universality class in two dimensions (not 
shown). Having performed the FSS analysis, these two critical points can be determined 
with a higher accuracy from the Binder cumulant crossing method^, as an intersection of 
the Binder parameter Uq curves for different lattice sizes L and O = Mgi {i = 1,2) (see 


Figs. 4(c) and 4(d)). The critical temperatures were determined as /cbTc/| J| = 0.7485±0.001 


at D/\J\ = 0 and fc^Tc/lJ| = 0.3663 ± 0.001 at D/\J\ = —2. Also the critical values of 
the Binder cumulants Uo{Tc) = 0.611— conhrm the Ising universality class. Furthermore, in 
Figs. 4(e)[4(f) we show that for the critical temperatures determined by the Binder cumulant 
crossing method the data for different lattice sizes indeed collapse on a single curve. The 
apparent hrst-order character of the phase transition at D/\J\ = —1.48 will be discussed 
below. 

Let us now examine the transition between the two low-temperature ferrimagnetic phases 
FR^ and FR^. Since the expected phase boundary is almost vertical to the x-axis, instead 
of the temperature dependencies of various thermodynamic functions it is more convenient 
to look into their single-ion parameter dependencies at a hxed temperature. By plotting 
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FIG. 6: (Color online) (a) Energy distributions for D/\J\ = —1.47 and different L. The respective 
temperatures are tuned by the reweighing technique to achieve approximately equal peak heights, 
(b) FSS analysis of the susceptibility xm^i the specific heat c, for D/\J\ = -1.47,-1.48 and 
— 1.4825. For D/\J\ = —1.48 and —1.4825 the log-log plots are respectively fitted to the exact values 
of the tricritical exponents ratios yt/t't = 1.85, at/t't = 1.60 and the exponent 2, corresponding 
to the system volume, (c) FSS analysis of DImsh D2ms\ for D/\J\ = —1.46, with the 

coefficients of determination E? for the respective fits (from top to bottom) 0.9999, 1.0000, 0.9999. 

the order parameters as increasing and decreasing functions of D/\J\ one can observe their 
discontinuous character and the appearance of hysteresis loops, the widths of which increase 
with decreasing temperature. This behavior is demonstrated in Fig. [5] for the order parame¬ 
ter rrtsi and L = 48. Such behavior signals hrst-order phase transitions. They seem to persist 
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FIG. 7: (Color online) Time evolution of the sublattice magnetizations at the points (a) 
(L>/|J|,A:bT/|J|) = (-1.487,0.33) and (b) (-1.47,0.35), for L = 48. 

even to higher temperatures at which the hysteretic behavior in not apparent any longer. 
The highest temperatures at which we still could observe some signs of first-order phase 
transitions, such as bimodal energy distribution, was ksT/\J\ = 0.35 and the corresponding 


value of D/\J\ = —1.47 (see Fig. 6 (a)) 


Nevertheless, the energy barrier separating the two peaks upon initial increase seems 
to decrease for larger L, which indicate that the phase transition may not be truly hrst 
order. Indeed, when we checked whether the specific heat and susceptibility scale with the 
system volume, as it should be in case of a first-order transition, we found that within 
the used lattice sizes the linear ansatz could not be established ioi D/\J\ = —1.47, as 


shown in Fig. 6 (b) We note that such a behavior that can lead to misinterpretation of a 
second-order transition as first order was also observed in the Blume-Capel model as well 
as in the frustrated Ji — J 2 Ising antiferromagnet on a square lattice^*^. Clearly first-order 
scaling is observed only at D/\J\ = —1.4825. Nevertheless, fairly good linear hts were also 
achieved in the case of D/\J\ = —1.48 with the exponents that are between 2 and the exact 
tricritical values (1.85 for 7 / 1 / and 1.6 for suggesting that the system is close to 


the tricritical point. On the other hand, for D/\J\ = —1.46, the FSS presented in Fig. 6 (c) 


clearly indicates a second-order phase transition. Thus, we roughly estimate the tricritical 
point at {Dt/\J\, fc^Ti/l J|) = (—1.47 ± 0.01, 0.35 ± 0.01). 

It is interesting to notice that the phase transition at this point is no longer between 


13 







































0.9 


0 . 8 - 
0.7- 
0 . 6 - 
^ 0.5 
0.-1 
0.3- 


.ee 


e®' 


)0ffi0O 


^ "’■'^3eeee§0^r-'’'’‘"P 


FR 


, IXI 
0.2- FR^ i>i<] 

> ' <1 

' |-<-first order 


-3 


-2 


-1 0 
D/|J| 


FIG. 8: (Color online) Phase diagram of the model I in (/c^T/l J| — D/\J\) parameter space. The 
empty circles represent the phase transition temperatures /csTc/l J| between the paramagnetic state 
P and the ferrimagnetic states FR^ (±1/2, ±1/2, =pl) and FR 2 (±1/2, ^1/2,0), estimated from the 
specific heat peaks for L = 48, the empty triangles mark the hysteresis widths at first-order 
transitions between the phases FR^ and FR 2 with the expected phase transition boundary marked 
by the dash-dot line. The filled symbols at finite temperatures show more precise values obtained 
from the FSS analysis and the Binder cumulant crossing, where the diamond is the tricritical point 
(TCP), the hexagon is the critical endpoint (CE) and the square at {Dc/\J\, /csTc/l J|) = (—3/2,0) 
represents the exact value of the GS transition point. 


the two ferrimagnetic phases FR^ and FR 2 bnt between the phase FR^ and the param¬ 
agnetic phase. This is in line with the above observation of the hrst-order-like featnres 
of the phase transition at D/\J\ = —1.4825 between the paramagnetic and FR^ phases. 
The fact that there are hrst-order phase transitions between FR^ and FR 2 as well as FR^ 
and paramagnetic phases can be verihed by looking at the relevant order parameters. For 
demonstration, in Figs. [7] we show segments of snblattice magnetization time series ob¬ 


tained at {D/\J\,kBT/\J\) = (—1.487,0.33) (Fig. 7(a)) and (—1.47,0.35) (Fig. 7(b)). In 
both cases we can see discontinnons switching between two phases, however, in either 
point those phases are different. Namely, at (—1.487,0.33) the snblattice magnetizations 
(mA, rriB, me) switch between the valnes characteristic for the states (±1/2, ±1/2,0) 
and (±1/2, ±1/2, ±1), while at (—1.47,0.35) they switch between the valnes characteris- 
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tic for the states (0,0,0) and (±1/2, ±1/2, =f 1). This hnding implies the existence of a 
critical endpoint (CE) at which the second-order transition bonndary between the para¬ 
magnetic and FR 2 phases joins the above discnssed hrst-order transition line at abont 
{D^J\JlkBTj\J\) = (-1.485 ± 0.005, 0.335 ± 0.005). 

The resnlting phase diagram is presented in Fig. [HI The empty circles represent the 
psendo-critical points determined from the specihc heat maxima at L = 48 and the empty 
triangles mark the metastable branches of the hrst-order phase transitions obtained from 
the order parameter hysteresis loops. The hlled circles and sqnares at hnite temperatnres 
represent respectively second- and hrst-order transition points obtained from the FSS anal¬ 
ysis and the hlled diamond marks approximate location of the tricritical point. The hlled 
sqnare at {Dcl\J\, /c^Tc/l J|) = (—3/2, 0) shows the exact locations of the gronnd-state phase 
transition. The expected hrst-order phase transition bonndary, marked by the dash-dot line, 
is obtained by a simple linear interpolation between the estimated tricritical and the exact 
GS transition points and only serves as a gnide to the eye. We note that in this highly snp- 
pressed mixed-phase region the hrst-order phase bonndaries can be located qnite precisely, 
for example, by mnlticanonical MC simnlations^*^. However, in onr case, we are only in¬ 
terested in approximate location of the phase bonndaries. Then, considering the exact valne 
of the GS transition point, the estimate location of the tricritical point and assnming no 
anomalons behavior, snch as reentrance, the low-temperatnre part of the phase bonndary 
mnst be practically vertical to the D/\J\ axis. 


2. Mode///; S = (1/2,1,1) 


Also for the model II the phase bonndary as a fnnction of the single-ion anisotropy 
parameter D/\J\ is estimated from the specihc heat maxima for L = 48, except for D/\J\ 
close to the critical valne of Dc/\J\ = where the bonndary becomes almost vertical. 

The only identihed LRO phase is the ferrimagnetic phase FR^^: (±1/2, ±1, =pl), present 
for D > He, and the phase transition is second order complying with the standard Ising 


nniversality class. This is illnstrated in Fig. 9(a) for a selected valne of H/| J| =0, in which 


we show that the obtained ratios of the critical exponents qM^s/^^Mss are in a 

good agreement with the 2D Ising nniversality class valnes 1 and 7/4, respectively. We also 
checked the consistency of the specihc heat critical exponent valne a = 0 by verifying the 
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FIG. 9; (Color online) (a) FSS analysis of the critical exponent ratios XjvM^^ and 'yUssl^Msi-: (b) 
the fourth-order cumulant Ums^ temperature dependencies for different L, and (c) the Ums^ data 
collapse analysis at D/|J| = 0. In (a) the coefficients of determination E? for the respective fits 
(from top to bottom) are 0.9998, 0.9999 and 0.9999. Inset in (b) shows an alternative way of the 
critical temperature estimation from the FSS analysis. 


logarithmic scaling (not shown). The critical temperatnre for D/\J\ = 0 estimated by the 


Binder cnmnlant method (Fig. 9(b)) and FSS analysis (inset in Fig. 9(b)) takes the valne 
/c^Tc/l J| = 1.0635 ±0.0015 and the cnmnlant cnrves for different L intersect at the nniversal 
valne of UM^siTc) = 0.611. 

On approach to the critical value Dc/\J\ = —3/2 the phase boundary rapidly drops and 
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FIG. 10: (Color online) Order parameter as a function of the increasing (>) and decreasing 
(<) single-ion anisotropy parameter D/\J\ at various temperatures and L = 48. 


becomes almost vertical. Therefore, in order to locate the critical temperatures in this re¬ 
gion, it is more convenient to measure the physical quantities at a fixed temperature as 
functions of the parameter D /\J\. At sufficiently low temperatures the measured quantities 
show some properties typical for hrst-order phase transitions. Namely, as the anisotropy 
parameter D/\J\ is decreased and increased at the fixed temperature the sublattice mag¬ 
netizations, the order parameter and the internal energy show discontinuities at some 
values of D /\J\^ as demonstrated in Fig. [10] for nisz- Nevertheless, it is interesting to notice 
that, in contrast to the strongly hysteretic behavior of the model I at the transition between 
the two ferrimagnetic phases, no apparent hysteresis can be observed in the present model. 
Discontinuous character of the transition reflected in bimodality of the relevant observables, 
such as the internal energy, disappears at higher temperatures but is still evident at temper¬ 


atures slightly above kBT/\J\ = 0.2. In Fig. 11(a) it is demonstrated for D/\J\ = —1.48 and 
the temperatures /cbT/|J| 0.215 tuned by the reweighing technique for each L to achieve 
distributions with the two modes of about the same heights. The inset gives us some idea 
about the characteristic tunneling times between the coexisting phases for L = 48. However, 
similar to the situation in the model I presented in Fig. [6l for D/\J\ = —1.48 the specific heat 
and staggered susceptibility maxima do not scale with volume. Such a scaling is observed 
only at slightly lower temperatures for D/\J\ = —1.4825 (at least for L > 72), as shown in 


Fig. 11(b) At higher temperatures the second-order phase transition with standard Ising 
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FIG. 11: (Color online) (a) Energy distributions for D/\J\ = —1.48 and different L. The respective 
temperatures are tuned by the reweighing technique to achieve approximately equal peak heights. 
Time evolution of the internal energy shown in the inset demonstrates tunneling between the 
coexisting phases for L = 48. (b) FSS analysis of the susceptibility and the specific heat 

c, for D/\J\ = -1.47,-1.48 and —1.4825. For D/\J\ = —1.4825 the log-log plots are fitted to 
the exponent 2, corresponding to the system volume, (c) FSS analysis of xm^sj D2ms3 

D/\J\ = —1.46 with the fitted Ising values of the critical exponents ratios = 1.75, l/vj = 1.00. 


critical exponents is recovered for D/\J\ = —1.46, althongh again larger system sizes are 


reqnired to reach the linear asymptotic regime (Fig. 11(c)). Thus, the tricritical point in 
the model II is roughly located at {Dt/\J\^ kBTt/\J\) = (—1.47 ± 0.01, 0.27 ± 0.04). 

The above results can be summarized into the phase diagram shown in Fig. [121 As in 
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FIG. 12: (Color online) Phase diagram of the model II in (/c^T/lJl — D/\J\) parameter space. The 
empty circles represent the phase transition temperatures kBTc/\J\ between the paramagnetic P 
and the ferrimagnetic phase (0, ±1,^1) estimated from the specific heat peaks for L = 48, 
the filled circle shows a more precise value obtained from the FSS analysis and Binder cumulant 
crossing, the filled diamond is the tricritical point and the filled squares represent a hrst-order 
transition point determined from FSS analysis at D/\J\ = —1.48 and the exact value of the GS 
transition point at Dc/\J\ = —2>l2. The empty triangles mark the first-order transition related 
discontinuities in the order-parameter ms 3 in the D/\J\ increasing (o) and decreasing (<i) processes. 

the phase diagram of the model I, the empty circles represent the pseudo-critical points 
determined from the specific heat maxima at L = 48 and the empty triangles mark the first- 
order phase transitions located from the jumps in the order parameter loops obtained by 
increasing and decreasing of the single-ion parameter D/\J\. The hlled circle at D/\J\ = 0, 
the square at Dl\J\ = —1.48 and the hlled diamond represent respectively second-order, 
hrst-order and tricritical points, obtained from the FSS analysis. As in Fig. |H1 the hlled 
square at {Dc/\J\,kBTc/\J\) = (—3/2,0) represents the exact value of the ground-state 
phase transition point. 


IV. CONCLUSIONS 

We have studied the mixed spin-1/2 and spin-1 Ising ferrimagnets on a triangular lattice 
with sublattices A, B and C, in two mixing modes: (Fa, Fb, Sq) = (1/2,1/2,1) (model I) and 
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(S'a, 5'b, Sc) = (1/2,1,1) (model II). The pure spin-1/2 and spin-1 Ising antiferromagnets on 
a triangular lattice show respectively no LRO^ and partial LRO for some range of a single¬ 
ion anisotropy parameter at low temperature with quasi-LRO of the Berezinskii-Kosterlitz- 
Thouless type at higher temperatures^. In comparison with these models in the present 
ferrimagnetic models the frustration is partially accommodated by different ferrimagnetic 
spin arrangements and thus their critical behavior is rather different from the pure systems. 
On the other hand, the net magnetization of both ferrimagnetic models is always zero and 
thus they can show no compensation points. Therefore, from this point of view, the behavior 
of the present mixed-spin models is typical for antiferromagnets rather than ferrimagnets. 

As for the critical properties, the model I shows two ferrimagnetic phases FR^: 
(±1/2, ±1/2, =pl) and FR^: (±1/2, =pl/2, 0), which may be compared with the ferromag¬ 
netic case displaying two ferromagnetic phases (±1/2, ±1/2, ±1) and (±1/2, ±1/2, 0}^. We 
note that on bipartite lattices the thermodynamic behavior of the systems with ferrimag¬ 
netic and ferromagnetic interactions is the same and thus the phase diagrams would be 
identical. However, there are substantial differences between the two phase diagrams for the 
triangular lattice ferromagnetic and ferrimagnetic models. First of all, due to frustration 
the transition temperatures for the ferrimagnetic case are signihcantly reduced and collapse 
with the ferromagnetic boundary only in the large negative D/\J\ limit, when the partial 
frustration inducing magnetic states on the C-sublattice are completely suppressed and the 
critical temperature for either case tends to the exact spin-1/2 Ising value on a honeycomb 
lattice ksTc/lJl = 0.3797—. The frustration is further increased close to the boundary be¬ 
tween the two ferrimagnetic phases FR^ and FR^, which is reflected in the depression in the 
order-disorder phase boundary that is absent in the ferromagnetic model. Nevertheless, the 
frustration did not seem to affect the standard Ising values of the critical exponents. We 
note that besides the above presented points of D/\J\ = —2 and 0, we also performed the 
FSS analysis in this region of an increased frustration at the P-FR^ branch of the phase 
diagram for D/\J\ = —1.6 (not shown) but did not hnd any deviation larger than statis¬ 
tical errors from the standard values. However, the most conspicuous difference from the 
ferromagnetic case is the presence of the strongly discontinuous phase transition between 
the ferrimagnetic phases (±1/2, ±1/2, =f 1) and (±1/2, ±1/2, 0), which is completely absent 
between the ferromagnetic phases (±1/2, ±1/2, ±1) and (±1/2, ±1/2, 0). This can be ex¬ 
plained by the fact that in the ferrimagnetic case dramatic changes occur when the two 
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spin-1/2 sublattices A and B, forming a connected honeycomb backbone, switch their mag¬ 
netizations between 1/2 and —1/2 and the spin-1 sublattice C between magnetic ±1 and 
nonmagnetic 0 states. On the other hand, in the ferromagnetic case nothing happens in 
sublattices A and B and in sublattice C the change between magnetic ±1 and nonmagnetic 
0 states occurs only gradually in the isolated (mutually directly noninteracting) spins. 

The model II has been shown to display only one ordered ferrimagnetic state FR^^: 
(0, ±1, =f 1) with no LRO on sublattice A and an antiferromagnetic LRO on the remaining 
sublattices B and C forming a honeycomb lattice. Similar phase occurs in the pure spin- 
1 triangular antiferromagnet for the single-ion anisotropy parameter —3/2 < D/\J\ < 0, 
however, it is destabilized for D/\J\ > 0^. On the other hand, in the present model II 
the increasing D/\J\ stabilizes the ferrimagnetic phase FR^^ and the critical temperature for 
D/\J\ —)■ oo tends to fc^Tc/l J| = 1.5188^, i.e., the exact value of the spin-1/2 Ising model 
on a honeycomb lattice when the spin states ±1 are considered in the Hamiltonian instead 
of ±1/2. In comparison with the ferromagnetic model II, again the critical temperatures 
are lowered due to frustration but qualitatively the phase diagrams look similar. As the 
parameter D /\J\ is decreased both models show change of the phase transition nature from 
second to hrst order at a tricritical point. Nevertheless, interestingly, the strong hysteretic 
behavior observed in the ferromagnetic model accompanying the first-order phase transitions 
is not evidenced in the ferrimagnetic one. 
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